      subroutine simple_dfirm1d(   n_fx, ny, n_gx, f, g, h)
      double precision f(0:(n_fx-1),ny)
      double precision g(0:(n_gx-1))
      double precision h(0:(n_fx+n_gx-2),ny)
      integer n_fx, n_gx, ny
c-----------------------------------------------------------------------------
c
c   This Fortran Subroutine written by
c	Jean-Pierre Panziera
c	Silicon Graphics Inc
c	September 12, 1991
c
c-----------------------------------------------------------------------------
	integer j, ihx, ifx, igx
	double precision zero
	parameter ( zero = 0.)
c-----------------------------------------------------------------------------
c
c   Compute the N 1D convolutions
c
c-----------------------------------------------------------------------------
	do j = 1, ny
	    do ihx = 0, n_fx+n_gx-2
		h(ihx,j) = zero
	    end do
	    do ifx = 0, n_fx-1
		do  igx =  0, n_gx-1
		    h(ifx+igx,j) = h(ifx+igx,j) + f(ifx,j) * g(igx)
		end do
	    end do
	end do

      return
      end


